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Interactions of different types of topological defects can play an important role in the aftermath 
of a phase transition. We study interactions of fundamental magnetic monopoles and stable domain 
walls in a Grand Unified theory in which SU(5) X Z 2 symmetry is spontaneously broken to SU(3) X 
SU(2) x U(l)/Z 6 . We find that there are only two distinct outcomes depending on the relative 
orientation of the monopole and the wall in internal space. In one case, the monopole passes through 
the wall, while in the other it unwinds on hitting the wall. 


I. INTRODUCTION 

Grand Unified Theories (GUTs) are based on large 
symmetry groups, the smallest of which is an SU(5) 
model with an additional, possibly approximate, Z2 sym¬ 
metry. When such large symmetries are broken in a cos¬ 
mological setting, several kinds of topological defects can 
be produced. The ensuing cosmology will depend criti¬ 
cally on the interactions of the different defects. In par¬ 
ticular, the SU(5) x Z2 symmetry breaking leads to the 
existence of magnetic monopoles and domain walls in the 
aftermath of the phase transition. We expect the mag¬ 
netic monopoles to interact with domain walls, poten¬ 
tially resolving the magnetic monopole over-abundance 
problem [1]. To investigate this idea further, we study the 
interactions of SU(5) monopoles and Z2 domain walls in 
this paper. 

The interaction of monopoles and domain walls was 
also studied in [2] with the domain wall structure given 
by 

$ = tanh (-) $ 0 (1) 

where the order parameter <£> is in the adjoint represen¬ 
tation of SU(5), is its constant vacuum expectation 
value (VEV), and w is the width of the domain wall. By 
numerical evaluation it was found that monopoles hit¬ 
ting this domain wall will unwind and spread on the wall. 
Subsequently, however, it was found [3-6] that the model 
actually has several domain wall solutions, including the 
one in Eq. (1), and that the lightest (stable) wall has a 
different structure (see Sec. II B). Hence the interaction 
of the stable wall and the monopole needs to be revisited. 

In Sec. II we provide details of the SU(5) x Z2 model, 
the monopole solution, the wall solutions, and finally our 
scheme for setting up a configuration with a monopole 
and a domain wall together. This provides us with initial 
conditions that we numerically evolve in Sec. III. The 
complexity of the field equations and the problem re¬ 
quires some special numerical techniques that we briefly 
describe in Sec. III. 

Our results are summarized in Sec. IV. Essentially we 
find that there are two internal space polarizations for 
the monopole with respect to the wall. One of the polar¬ 
izations is able to pass through the wall with only some 


kinematic changes. The monopole with the other polar¬ 
ization is unable to pass through the domain wall and 
unwinds on the wall, radiating away its gauge fields. The 
disappearance of this monopole is further explained in 
Sec. IV. 


II. THE MODEL 


The SU(5) model we consider is given by the La- 
grangian: 

L = Ax a ^X a ^+^D lx (j) a D tl (t> a -V(^) ( 2 ) 

where <£> = (f a T a (a = 1,..., 24), X“ u are the gauge field 

strengths defined as 

= d^X v - dyX^ - ig[X„, X v \ , (3) 

Xn = X^T a are the gauge fields and g is the coupling 
constant. T a are the generators of SU(5) normalized by 
Tr(T a T 6 ) = 5^/ 2. The covariant derivative is given by 

D^=d^ a -ig[XnM a - ( 4 ) 

The most general renormalizable SU(5) potential is 
V(ff>) = —m 2 Trff> 2 +7Trff> 3 +h(Trff> 2 ) 2 +ATrff >4 —Vo , (5) 


and we will assume that 7 vanishes, giving the model an 
additional Z 2 symmetry. For A > 0 and h + 7A/30 > 0, 
the potential has its global minimum at [7] 

<*>o = -^=diag(2,2,2,-3,-3), (6) 


with r\ — mf yjh + 7A/30. The VEV, 4>o, spontaneously 
breaks the SU(5) symmetry to SU(3) x SU(2) x U(l)/Z 6 . 

In what follows, the four diagonal generators of SU(5) 
are chosen to be 


A 3 = ^diag(l, —1, 0,0, 0), 

A 8 = ^=diag(l, 1,-2,0,0), 
r 3 = ^diag(0, 0,0,1,-1), 
y=2=diag(2,2,2,-3,-3). 


( 7 ) 
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We use a = 1,2,3 to denote generators T a = r a = 
diag( 0 , 0 , 0 , cr a / 2 ) where a a are the Pauli spin matrices. 

A. The monopole 

Let us consider a magnetic monopole whose winding 
lies in the 4-5 block of 4>. This is possible [ 8 ] if we take 
the VEV along one of the radial directions far away from 
the monopole to be 

$°° = ^-^|^diag(2, —3, 2, 2, —3) 

= V^(h+ T 3) + l( Y -V5X 8 ). (8) 

The monopole ansatz for the scalar field can be written 
as [3] 

$M(r) = P(r) x a Ta+M(r) ^A 3 - ^A 8 j +N(r)Y, 

( 9 ) 

while the non-zero gauge fields can be written as 

X i=^^~K(r)),(a = l, 2,3) (10) 

and P(r), M(r), N(r), and K(r ) are profile functions 
that depend only on the spherical radial coordinate r = 
x 2 + y 2 + z 2 and satisfy the boundary conditions: 

/ 5 y/b 

—, M(oo) = T)—, 

1Z o 

N(oo) = K(oo) = 0. (11) 

The profile functions for the monopole alone were evalu¬ 
ated numerically and are shown in Fig. 1. 

The non-Abelian magnetic field can be defined as [9] 

B k = 

with the associated energy density given by Tr(BkB k ). 
Far away from the centre, the monopole field becomes 
B k Qx k /(gr 3 ), with Q = r^x^ jr. 

The monopole charge Q includes a component along 
the generator of the unbroken U(l) symmetry (4>oo of 
Eq. ( 8 )), as well as SU( 2 ) and SU(3) magnetic charges. 
The U(l) part of the magnetic field, which is a defining 
feature of a topological SU(5) monopole, is given by 

By = -~e ijk X^ a ( 12 ) 

2 J 

where = (j) a /yf (j) b (j) b . As discussed in [10], other def¬ 
initions of the Abelian magnetic field are possible, and 
these differ from our definition but only within the core 
of the monopole. Since we only use our definition to plot 
the long range Abelian magnetic field (see Fig. 3) the 
definition in Eq. (12) is sufficient. 



FIG. 1. The profile functions for the monopole alone, evalu¬ 
ated numerically, for a model with r) = 1, h/X = —0.2 and 
A = 0.5. 

B. The wall 

Without loss of generality [4], the domain wall solution 
can be taken to be diagonal at all 2 ) and written in terms 
of the diagonal generators of SU(5) as 

®dw{z) = a(z)X 3 + b(z) A 8 + c(z)r 3 + d(z)Y . (13) 

In each of the two disconnected parts of the vacuum man¬ 
ifold M there are a total of 10 different diagonal VEVs 
corresponding to all possible permutations of 2’s and 3’s 
in Eq. ( 6 ). Topology dictates that there must be a do¬ 
main wall separating any pair of VEVs from the two dis¬ 
connected parts of M. However, not every such pair of 
VEVs corresponds to a stable domain wall solution. For 
instance, as shown in [3], the wall across which 4>o goes 
to — 4>o is unstable and will decay into a lower energy 
stable wall. The stable domain walls are obtained when 
both 3’s in Eq. ( 6 ) change into 2’s across the wall. 

Let us choose the boundary condition at z = — oo to 
be 

= $(z = -oo) = ^=diag(2, -3, 2, 2, -3) 

= + T 3)+ ^0^ ~ V5Xs). (14) 

For this choice of <£>_, there are three different choices of 
4>(z = Too), proportional to 

diag(3, —2, —2, 3, —2) 
diag(—2, -2, 3, 3,-2) 
diag(3, —2,3, —2, —2), (15) 

that lead to stable domain walls. For the purpose of un¬ 
derstanding the monopole-wall interactions, it is suffi¬ 
cient to consider only two of the above, corresponding to 
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Case 1 - Case 2 



FIG. 2. The wall profile functions for Cases 1 and 2 for a 
model with 77 = 1, h/X = —0.2 and A — 0 . 5 . Note that the 
profile function c(z) goes to zero in Case 2, which gives an 
unbroken SU( 2 ) C SU( 3 ) symmetry in the 4-5 block of 
The profile function a(z) is the same for Cases 1 and 2. 


that the symmetry in the 4-5 block is different for the 
( 2 ) 

and vacua. This is going to be of direct relevance for 
the fate of the monopoles. 

The profile functions a(z), b(z ), c(z) and d(z) for both 
cases are shown in Fig. 2. In each case, they are linear 
combinations of two functions F+(z) and F_(z) defined 
by the alternative way of writing the domain wall solution 

[3] 


_ *+(*) - M*) T? M , M*) + M*) 7? , X 

®DW ~ -o- *~\ z ) H-o- * + \ z ) > 


( 18 ) 

where F + (±oc) = 1, F_(±oc) = =tl. For a general choice 
of parameters, functions F±(z) must be found numeri¬ 
cally. For h/X = —3/20, they are known in closed form 
[3]: F+(z) = 1, F-(z) = tanh(m^/v // 2 ). Correspondingly, 
for this value of h/X, the four functions a(z), b(z ), c(z), 
and d(z) are either constant or describe a transition from 
one constant value to another. For hj A / —3/20 the 
“constant” functions develop a small bump around z = 0 
as can be seen in Fig. 2. 


C. Monopole and Wall 


the two distinct entries in the 4-5 block of <F. We take 
the first to be the same as in [3], subsequently referred 
to as Case 1 : 

A 1} = ^7ff diag ( 3 ’ _2 ’ _2 ’ 3 ’" 2 ) 

= + r 3 ) - \/5A 8 ). (16) 

The value of the field in the core of this wall is pro¬ 
portional to diag(l, —1, 0,1, —1). The other case, subse¬ 
quently referred to as Case 2 , has 


As our initial configuration, we take the monopole to 
be on the z = —oo side, far away from the wall. In this 
case, the ansatz for the initial combined field configura¬ 
tion of the wall and the monopole can be written as [3] 


Sm+dw = P{r) N x a Ta + N(v)^\y 
c(— oo) di— oo) 

v 7 a= 1 v 7 


+ M(r) 


\/3 a(z') 


A.3 - 


2 a(— oo) 2 b(— oo) 


(19) 


$ + 2) = 27n dias(3, _2 ’ 3, _2, _2) 

= + |(4F - v^As), (17) 

with the field in the wall being proportional to 
diag(l, —1,1,0,—1). A novel feature of these walls is that 
the unbroken symmetry groups on either side of the wall 
are isomorphic to each other but they are realized along 
different directions of the initial SU(5) symmetry group. 
Hence the wall is the location of a clash of symmetries 
[ 12 ]- 

Note that the symmetry within the wall is [SU(2) x 
U(l)] 2 . The SU( 2 )’s correspond to rotations in the 1-3 
and 2-5 blocks and the U(l)’s to rotations along as in 
the 1-2 and 3-5 blocks. Therefore the symmetry group 
within the wall is 8 -dimensional, and is smaller than the 
12 -dimensional symmetry outside the wall 1 . Also note 


where z' = ^{z — Zo), 7 = 1/y/l — v 2 is the boost fac¬ 
tor, v is the wall velocity and zq is the initial posi¬ 
tion of the wall. The monopole isatx = 0 = y = z. 
It is easy to check that, far away from the monopole, 
the profile functions take on the values in Eq. (11) and 
$m+m ^dw- Close to the monopole, z r -T — 00 , 
since the monopole is initially very far from the wall, 
and —>■ as desired. We work in the tempo¬ 

ral gauge, Xg = 0 , and with the initial ansatz for the 
gauge fields given by Eq. (10) for both cases. 

It is instructive to examine the difference in the nature 
of the magnetic field in Cases 1 and 2 . As mentioned 
in Sec. II A, the charge of our monopole along the z- 
direction, Q = (l/2)diag(0, 0, 0,1, —1), is a combination 
of the U(l), the SU( 2 ), and the SU(3) magnetic charges. 
Since the VEV of <f> in our model is along the generator 


1 For simplest domain walls, such as kinks in A<f> 4 , the full symme¬ 


try of the Lagrangian is restored inside the core. However, the 
symmetry inside stable domain walls in SU(IV) X Z 2 is always 
lower than that of the vacuum [4]. 
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of (hypercharge) U(l), the magnetic field, as defined in 
Eq. (12), corresponds solely to the U(l) component of the 
charge. In Case 1, Tr(Q4>) is the same on both sides of 
the wall and the U(l) magnetic field is unaffected by the 
presence of the DW. In Case 2, however, Tr(Q4>y ; ) = 
0 and there is no magnetic field corresponding to the 
unbroken U(l) on the z = Too side of the wall. Instead, 
the gauge field on that side is associated with an SU(2) 
subgroup of the unbroken SU(3). We note that, while the 
magnetic energy density associated with the gauge field 
is unaffected by the presence of the wall, it is specifically 
the U(l) magnetic field that is a defining feature of a 
topologically stable monopole. 


Case 1 Case 2 
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FIG. 3. The magnetic field B y (defined in Eq. (12)) multiplied 
by r 2 for Cases 1 and 2, where at each point r 2 By and r 2 By 
are plotted as a vector. In Case 1, there is a magnetic field 
associated with the unbroken U(l) symmetry on both sides of 
the wall. In Case 2, the magnetic field becomes associated with 
the SU(2) C SU(3) on the z = +oo side on the wall, while its 
U(l) component vanishes. Note that it is the U(l) magnetic 
field that characterizes a topologically stable monopole. 

The magnetic field, as defined in Eq. (12), is plotted for 
both cases in Fig. 3, where the vectors have components 
r 2 By and r 2 By • This plot shows that, in Case 1, there 
is a U(l) magnetic field on both sides of the wall falling 
off as r 2 as expected, while in Case 2 the U(l) magnetic 
field is zero on the z = Too side of the wall. 


III. EVOLUTION 

Let us consider an initial monopole-wall configuration 
given by Eq. (19) in which VEV at z = — oc is given by 
4>_ in Eq. (14). As mentioned in the previous Section, 
there are 2 types of boundary conditions at z = Too, 
given by Eqs. (16) and (17), dubbed Case 1 and Case 
2, leading to 2 different outcomes of the monopole-wall 
collision. 

Before considering the two cases in detail, let us note 
that initially, when the monopole and the wall are very 
far away from each other, the field configuration has just 
three non-zero gauge fields and six scalar fields corre¬ 
sponding to the generators that appear in Eq. (19). Be¬ 


cause these six generators form a closed algebra, it follows 
from the equations of motion that the subsequent evolu¬ 
tion does not involve fields corresponding to the other 18 
generators. Namely, the scalar and the gauge field equa¬ 
tions are 

D^(j) a = -dV/d(j) a (20) 

D^X^ a = g f abc (D v <&) b <j) c (21) 

where f abc are the SU(5) structure constants defined by 
[T a ,T b ] = ifabcT c . Let C be the set of indices of the 6 
generators that appear in the initial field configuration 
given by Eq. (19). Since the 6 generators form a closed 
algebra, f abc = 0 for a £ C and 6, c G C. Now let <fi a and 
X® be fields corresponding to any a £ C. If </> a and X are 
zero at the initial time, they will remain zero if f abc = 0 
for 6, c G C and dVjd(\) a ^ 0. The former condition is 
satisfied as mentioned above, while the latter holds since 
Tr [T a T b ] oc S ab and Tr [T a T b T c T d ] = 0 for 6, c, d G C, 
as we have checked by explicit evaluation. Thus, for our 
purposes, it is sufficient 2 to consider only a G C. 

Our numerical implementation is based on techniques 
developed in [2]. First, the DW and the monopole pro¬ 
file functions are found via numerical relaxation. The 
monopole is initially located at the center of the lattice. 
We give the DW a velocity towards the monopole and 
boosted profiles are inserted into the initial configuration 
given by Eq. (19). With the initial time derivatives sim¬ 
ply determined from the Lorentz boost factor, this initial 
configuration is evolved forward in time using a staggered 
leapfrog code. The boundary conditions require special 
care since the wall extends all the way across the lattice. 
We have implemented boundary conditions in which the 
field is extrapolated across the boundary. We have nu¬ 
merically tested that this boundary condition leads to 
a smoothly evolving domain wall, without any spurious 
incoming radiation. Even though our problem has axial 
symmetry, we work in Cartesian coordinates as this offers 
superior stability. However, as discussed in [13], we take 
advantage of the axial symmetry of our configuration to 
restrict the lattice to just three lattice spacings along the 
y direction. We then use a 256 x 256 lattice grid for the 
x and z coordinates. Additionally, the axial symmetry 
allows us to solve only for positive x and use reflection to 
find the fields at negative x. The units of length are set 
by r] = 1 and we take each lattice spacing to correspond 
to half of a length unit. In these units, the range of x and 
z axis for a 256 x 256 grid is [—64, 64]. Note that in some 
figures we do not plot the entire lattice. The radius of the 
monopole core is about 10 length units and is about the 
same as a half of the domain wall width. At the initial 
time, the wall is 30 length units away from the center of 
the monopole. 


2 Although the field components for a 0 C continue to vanish dur¬ 
ing evolution if they vanish initially, we cannot exclude the possi¬ 
bility that the fields in these other directions may grow unstably 
if they did not vanish initially. 
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Initial 


Final 



FIG. 4. The potential and magnetic energy densities in the xz plane for the colliding monopole and wall in Case 1. We see 
that the monopole passes through the wall and the energy densities remain localized. Additionally, we see the magnetic energy 
density is unchanged before and after the collision. 
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FIG. 5. The scalar field <f> a in the xz plane for the colliding 
monopole and wall in Case 1, where at each point (j) 3 and 
0 1 are plotted as a vector. In this case the scalar field ar¬ 
rangement in direction and magnitude remains virtually un¬ 
changed. 


A. Case 1: the monopole passes through 


It is not difficult to predict that the monopole in Case 
1 will pass through the wall. The monopole winding is 
due to the fields in the SU(2) subgroup corresponding 
to generators r a , a = In Eq. (19), these fields 

are multiplied by the function c(z ) which has the same 
value at z = Too and, as known from [3], is approxi¬ 
mately constant across the domain wall. Only b(z) and 
d(z) change signs across the wall, but these are irrelevant 
for the winding of the monopole. Thus, the presence of 
the wall is of no qualitative consequence to the wind¬ 
ing of the monopole or its profile functions. The only 
effect is the small change in c(z) around z = 0 (note 


that, as mentioned earlier, c(z) is strictly a constant when 
h/X = —3/20). 

We numerically collide the monopole and the wall by 
giving the wall an initial velocity of 0.8 (in speed of light 
units) and choosing parameters p = 1, h = —A/5, and 
A = 0.5 for F(<F). 

Fig. 4 shows the potential and magnetic energy densi¬ 
ties as the wall hits the monopole in Case 1. In addition, 
we plot the scalar field configuration in Fig. 5, where each 
point is a vector with components </> 3 and 0 1 . These fig¬ 
ures show that the magnetic energy density and the scalar 
field configuration remain unchanged after the collision, 
and that the potential energy densities corresponding to 
the monopole and the domain wall remain localized. This 
does not imply a complete absence of interaction between 
the wall and the monopole - some interaction is expected 
due the non-linearity of the scalar field potential. 

To see if the monopole gains momentum due to the in¬ 
teraction, we have evaluated the centre of energy (COE) 
defined as 


^coe(^) = 


f y d 3 x z p(£, x) 
fv d 3 x p(t, x) 


( 22 ) 


where V is the volume of a finite cylindrical region cen¬ 
tred at the origin and extending l/8th of the lattice size 
in the x- and ^-directions, while p is the energy density. 
For h/X = 0, we give the wall a velocity of v = 0.9 
towards the monopole and compare the initial zqeo to 
the one after the wall passes away. We see a very slow 
drift of the COE in the direction of the wall velocity. 
We performed the same procedure using different model 
parameters and wall velocities and the outcome was qual¬ 
itatively the same. In all cases, while the direction of the 
drift is clear, the magnitude is extremely small and too 
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close to the numerical uncertainties to allow a definitive 
quantitative analysis. 

B. Case 2: the monopole unwinds 

As in Case 1, it is possible to guess the outcome of the 
monopole-wall collision without doing numerical simula¬ 
tions. For this, we note that has an SU(2) symmetry 
in the 4-5 block, which means that there is no topology 
that can support the winding. Thus, the monopole cannot 
exist in that corner of the matrix. An equivalent way to 
see this is to note that the function c(z), which multiplies 
the three relevant monopole scalar fields, goes to zero at 
z = Too (see Fig. 2), effectively erasing the monopole. 

Additional insight can be gained by noting that the 
long range magnetic field of the monopole transforms into 
an SU(3) magnetic field on the far side of the wall. More 
explicitly, the U(l) magnetic field is given by Eq. (12) 
with Xfj determined using the solution in Eq. (10). Since 
Xij only has components in the r a directions, it lies in the 
4-5 block. However, the 4-5 block is entirely within the 
unbroken SU(3) on the right-hand side of the wall. Thus 
the long range magnetic field of the monopole is purely 
SU(3) on the right-hand side of the wall and, from the 
vantage point of someone there, there is no U(l) mag¬ 
netic field emerging from the left-hand side of the wall. 
However, a U(l) magnetic field is an essential feature of 
a topological monopole. Thus, from the right-hand side 
of the wall, there is no magnetic monopole in the system, 
only some source of SU(3) magnetic flux. 

Doing the numerical simulation with the parameters 
chosen as before, we plot the potential and magnetic en¬ 
ergy densities as the wall hits the monopole in Fig. 6. This 
figure shows that the potential energy for the monopole 
disappears as the wall and monopole collide, and the 
magnetic energy that was stored in the monopole ra¬ 
diates away in a hemispherical wave. The collision was 
simulated with initial wall velocities ranging from 0.1 to 
0.99 for h/X = —1/5, and initial wall velocities of 0.6, 0.8 
and 0.99 for h/X = —3/20 and 1/5. In all of these cases, 
the result of the collision was unchanged. 

In Fig. 7, we show the a = 1,2,3 components of the 
scalar field using two different representations. In the first 
row, the fields 0 3 and 0 1 are plotted as a vector. The plot 
shows that the components of the field that are responsi¬ 
ble for the winding vanish on the z = +oo side of the wall. 
In the second row of Fig. 7, the color represents the mag¬ 
nitude |0| = V0 a 0 a , a = 1,2,3, while vectors are drawn 
of fixed length and direction given by tan _1 (0 3 /0 1 ). Even 
though |0| becomes very small, it is not strictly zero at a 
finite distance from the wall, and so one can still define 
the direction of the arrow in this way. One can see that 
initially the field has a hedgehog configuration across the 
wall. However, as the wall sweeps along, the fields on the 
z = Too side of the wall rotate around in such a way as 
to unwind the monopole. In the final step, all fields that 
are non-zero are pointing in one direction, and therefore 


the monopole winding is gone. 

IV. CONCLUSIONS 

In a Grand Unified model there can be several types of 
defects, including magnetic monopoles and domain walls. 
In the aftermath of the cosmological phase transition in 
which the Grand Unified symmetry is spontaneously bro¬ 
ken to the standard model symmetry, the monopoles and 
walls will interact 3 . We have studied these interactions 
explicitly in an SU(5) x Z 2 GUT, taking into account that 
the model has several different types of domain walls, 
and that only the lowest energy wall is expected to be 
cosmologically relevant. Even this stable wall has several 
different orientations in internal space, two of which are 
distinct for the purposes of monopole-wall interaction. 

The first wall (Case 1 above) is found to be trans¬ 
parent to the monopole. This is simply because the do¬ 
main wall mainly resides in a certain block of field space, 
while the winding of the monopole resides in a differ¬ 
ent non-overlapping block. The interactions between the 
monopole and the wall are very weak, and only affect the 
dynamics of the monopole as it passes through the wall. 
Depending on the parameters, the monopole might be 
attracted or repelled by the wall leading to a time delay 
or advance as the monopole goes through. 

The second wall (Case 2 above) is opaque to the 
monopole. When the monopole hits the wall its energy is 
transformed into radiation on the other side of the wall, 
as seen in Fig. 6. A useful way to picture this system 
is to consider a magnetic monopole that is located in¬ 
side a spherical domain wall. Now there is a topological 
magnetic monopole inside the wall, but only an SU(3) 
magnetic flux from the outside. In particular, there is no 
topological magnetic monopole as seen from the outside. 
Therefore the spherical wall itself must carry the topo¬ 
logical charge of an antimonopole 4 . If the spherical wall 
shrinks, either it can annihilate the magnetic monopole 
within it and radiate away the energy, or the monopole 
can escape the wall, in which case the wall would then 
collapse into an antimonopole so that the total topolog¬ 
ical charge of the system continues to vanish. Our ex¬ 
plicit numerical evolution shows that annihilation occurs 
for the parameter ranges we have considered. We note 
that the unwinding of the monopole in the Case 2 may 
be related to the mechanism of formation of non-Abelian 
clouds (massless monopoles) [15]. 

Our results have bearing on cosmology as they explic¬ 
itly show the possible destruction of magnetic monopoles. 
In the case where the Z 2 symmetry is approximate, the 
walls will eventually decay away, and it is possible that 
these interactions could lead to a universe that is free of 


3 Scattering of fermions and GUT domain walls was studied in [14] 

4 The correspondence between spherical domain walls and global 
monopoles in SU(JV) has previously been noted in [5]. 
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FIG. 6. The potential and magnetic energy densities in the xz plane for the colliding monopole and wall in Case 2. We can see 
that as the domain wall and monopole collide, the potential energy contained by the monopole disappears and the monopole 
begins to radiate away its magnetic energy in a hemispherical wave. Note that the middle and final plots for the magnetic 
energy density have a much smaller scale as the ripples are not visible at the original scale. 


magnetic monopoles. Estimates in [1] indicate that this 
possibility is worth investigating in more detail. With 
several types of domain walls and monopoles simultane¬ 
ously forming in a phase transition [16-18], and with the 
complex nature of both the inter-wall [19] and monopole- 
wall interaction, the fate of the monopoles will remain 
uncertain until a comprehensive simulation of the GUT 
phase transition is performed. We leave this for a future 
study. 
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FIG. 7. The scalar field cj) a in the xz plane for the colliding monopole and wall in Case 2. At each point in the first row, (j) 3 
and 0 1 are plotted as a vector. In the second row, the length of the arrow is fixed, while the direction of the arrow is given by 
tan _1 (0 3 /0 1 ) and the color represents the magnitude of the field \<f>\ = V <j) a <j) a for a — 1, 2, 3. The first row shows the monopole 
unwinding as the wall sweeps past it, and the second shows how the fields arrange themselves to unwind the monopole. 

























